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TECHNICAL NOTE 2116 


LINEARIZED SUPERSONIC AXIALLY SYMMETRIC 
FLOW ABOUT OPEN -NOSED BODIES OBTAINED 
BY USE OF STREAM FUNCTION 
By Franklin Moore 


-SUMMARY 

The use of Stokes' stream function for axially symmetric prob- 
lems in linearized supersonic flow is discussed. The computation 
of streamlines away from the body is shown to be facilitated by 
solving the stream-f unction problem rather than the potential prob- 
lem. Use is made of the stream function throughout this report. 

Half -power supersonic source distributions along the axis of 
symmetry are shown to provide streamlines that, at finite distances 
from the axis, turn a sharp comer at the Mach cone and are regular 
downstream of the Mach cone. Distributions of this type are used 
to provide the following procedures for open-nosed bodies of revo- 
lution at zero lift; 

1. A numerical method for the computation of pressures along 
the body: The method, for which a simple example is given, is 
especially easy to apply on the forward part of the body. 

2. Linearized formulas for the variation of strength and angle 
of the attached bow shock wave as functions of radial distance 
from the axis . 


INTRODUCTION 

In the design of a supersonic aircraft, the aerodynamics of 
an open-nosed component (for example, a ram-jet -engine housing) 
must often be described. In the analytical treatment of this type 
of problem, the determination of the flow about slender, open-nosed 
bodies of revolution at zero yaw may frequently be of value. 

A standard method for computing the external pressure distri- 
bution associated with this type of flow is the numerical procedure 
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described in reference 1, wherein a supersonic source distribution 
on the axiff of symmetry is used to distort the free stream in such 
a way that one of the resulting stream surfaces approximately repre- 
sents the body surface. A source distribution automatically meets 
two requirements of the problem: 

1. The flow produced satisfies the linearized equations of 
motion of a compressible fluid. 

2 . Flow disturbances are limited to the region bounded by the 
outer surface of the body and the leading Mach wave proceeding out- 
ward and downstream from the lip, provided the source distribution 
is considered zero upstream of the apex of the Mach cone of which 
the leading Mach wave is a frustum. 

A particular source distribution must be selected to provide 
a stream surface coincident with the body, to the first approxi- 
mation in body slope. In the case of slender pointed bodies of 
revolution (reference 2), a simple relation exists between the 
source strength distribution and the slope of the body to be repre- 
sented. The existence of this relation is due to the nearness of 
the body surface to the axis; the streamline slope at a point near 
the axis is a consequence only of the local value of the source- 
distribution function, to the first approximation. 

When the body has an open nose, no such simplification may be 
considered to apply to its forward portion. The source distribution 
is then determined by an integral equation, which must be numeri- 
cally solved, inasmuch as no closed solution seems available. 

If the boundary condition on streamline slope at the body sur- 
face is applied on a mean cylinder, this integral equation may, in 
principle, be solved by operational methods, yielding solutions 
equivalent to those obtained for "quasi -cylindrical ducts" by 
Lighthill, whose work is summarized and extended in reference 3. 

The solutions appearing in reference 3 are given in integral form, 
and must be evaluated by numerical quadrature at each point where 
the pressure is to be determined. The accuracy of these results 
will suffer, in certain cases, from the application of boundary 
conditions at a mean cylinder. 

A form for the source distribution is assumed in reference 1 
(and in the present report) and the und etermined parameters are 
then selected in such, a way that boundary conditions are satisfied 
with adequate accuracy. The form assumed in reference 1 is trape- 
zoidal; that is, has a step-wise constant slope. The number of 
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steps and their spacing is determined "by the accuracy desired, and 
the slope of the distribution function between stations is deter- 
mined by the shape of the body to be represented. 

In selecting a form of solution wherein only certain param- 
eters are available for variation, it is desirable to anticipate, 
insofar as possible, the features of the "exact" solution. On 
this basis, the trapezoidal distribution appears to have a certain 
disadvantage: namely, when approaching the intersections of the 
body surface with the Mach cones emanating from the points on the 
axis where the source-distribution slope is discontinuous, the 
streamline curvature tends toward infinity, whereas the streamline 
slope remains finite and continuous and, in fact, vanishes at the 
leading Mach wave . Yerif ication of this statement is subsequently 
provided. This singular behavior of the streamline chosen to- 
represent the body' means that a large number of stations is 
required in regions of rapidly varying body slope, for example, 
near the lip of an open-nosed body. In particular, adequate evalu- 
ation of the local effect of discontinuous body slopes required a 
relatively dense local concentration of integration stations. An 
analysis of the linearized supersonic flow- over the forward part 
of such bodies was made at the HACA Lewis laboratory and is pre- 
sented herein. 

A family of solutions is developed for axially symmetric open- 
nosed bodies at zero yaw, which may be combined to provide either 
smooth or polygonal streamlines, as desired, and whose parameters 
can conveniently be selected in order best to represent the forward 
part of a given body. In particular, solutions are sought that can 
be smoothly matched to the body in the same manner as a polynomial 
could be. The supersonic nature of the flow allows the superposition 
of smooth solutions in such a way as to provide streamline corners. 

In brief, bodies encountered in practice will be smooth 
between and approaching corners . Solutions whose streamlines 
exhibit regular behavior near a leading Mach wave, but which may 
have corners precisely at the Mach wave, could therefore be matched 
to the forward portion of a given body with a minimum of "forcing" 
and hence a minimum of computational labor. 

TTse is made in this analysis of Stokes • stream function because 
solutions in terms of a stream function can readily be solved for 
streamlines, which may directly be compared with the given body 
contour. The stream function for linearized axially symmetric 
supersonic flow was used by Ferrari (reference 4), who developed 
the differential equation and solved it for pointed bodies of 
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revolution "by a Fourier integral method. In idle analysis that 
follows, the use of the linearized stream function is discussed in 
somewhat more detail than appears in reference 4. In the interest 
of continuity, certain of Ferrari' s results are rederived. The 
stream function is, of course, not essential to the analysis; 
equivalent results could he obtained by using the velocity 
potential . 


SYMBOLS 

The following symbols are used in this report: 


pressure coefficient, 


P-P_ 


p P u 
o oo 00 


E 

f(x) 

*n(n) 

K 

*n 

M 

N 

P 

Qn(n) 

B 


L E 


complete elliptic integral of second hind 

local strength of supersonic source distribution 

dimensionless functions corresponding to ¥, ¥^, and 
¥ aj , respectively 

flow function 

complete elliptic Integral of first kind 
constant' 

Mach number 

number of terms in terminated series 
static pressure 
flow function . 

radial coordinate of lip of open-nosed body 

radial coordinate measured from, and normal, to, axis of 
symmetry 

radial coordinate with which given streamline enters 
leading Mach cone 


dimensionless radial coordinate 
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t 

u 

Ufp 

U-oo 

T 

t t 

X 


7 

a 

P 

7 

5 

e(t) 

n 

Q, Sjj 

n 

V 

e • 

p 

ct(tj) 

T 


flaw function 

dimensionless axial coordinate 
axial perturbation velocity 
total axial -velocity component 
velocity of free stream 

» 

radial perturbation velocity 

total radial-velocity component 

axial coordinate measured in direction of axis of 
me try, origin at apex of leading Mach, cone 

variable of integration 

initial body slope 

cotangent of Mach angle of free stream, 

ratio of specific heats 

small quantity depending on body shape, 

body fineness parameter, Sg-1 

ratio of s to t 

angle between bow shock wave and axis of symmetry 
refers to dimensionless coordinates 

measure of distance from lip, x-pr-g 

measure of distance from lip, 1 -tj 

variable of integration 

density 

argument of certain elliptic integrals 
variable of integration 


rg-E 
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cp variable of integration 

¥(x,r) perturbation stream function 

¥ T (x,r) Stokes' stream function for supersonic axially symmetric 
flow 

¥ free -stream value of stream function 

00 

Subscripts : ’ 

0 stagnation value 

00 free -stream value 

B value on body 

1 index identifying point on body 
The prime denotes ordinary differentiation. 


ANALYSIS 

Stokes' Stream Function for Linearized 


Supersonic Axially Symmetric Flow 


Stokes’ stream function may be defined in the following manner 
(reference 5, par. 94) : 


Urp == 


1 ^T (x,r) 
P r dr 


r T = 


_ f3» I (x,r) 
P r Bx 


( 1 ) 


thus identically satisfying the exact equation of continuity 
(restricted to steady flow in this compressible case), which may 
be written 


s (o V + I k (rf V 


0 
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By comparison with the treatment of the incompressible case in 
reference 5, it is clear that, here also, d¥ T is an exact dif- 
ferential and represents the mass flux across a line element. 

Thus, a surface in space on which ¥rp is constant is a stream sur- 
face of the flow. 


In accordance with the assumptions of linearized theory. 



*T = ¥ » + ¥ 


( 2 ) 


where u,v « Uoo, and ¥<<¥«,. The introduction of these assump- 
tions into the isentropic gas law 



and the energy equation, written in the form 


p *0 

P Pq 



yields the approximate relation 


■f * 1 + M » 


2 u 


U, 


oo 


(3) 


Substituting equations (2) and (3) into definitions (1) and col- 
lecting terms of like order yields 


t d¥ 

1 CO 

03 r Sr 


u = - 


or ¥ ® = \ ^ r2 
1 S¥ 

p 1 2 r 


( 4 ) 

(5) 


1 S¥ 

‘ r 


( 6 ) 
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Differential equation and, 'boundary conditions . - Substitution 
of equations (5) and (6) into tbe irrotationality "" condition , 

5u 5 t . 

Sr Sz 


gives tbe following differential equation: 


«2 tl - r A- fl 

£x 2 Sr \r dr/ 


0 


( 7 ) 


Tbe flow over a given slender obstacle is obtained by solving dif- 
ferential equation (7), subject to tbe following boundary conditions: 


1. As in tbe case of tbe velocity potential, either a bound- 
ary condition at the bow Mach wave must be introduced or disturb- 
ances must be required to propagate only within downstream Mach 
cones, which is accomplished by use of the "supersonic source.". 

In order that streamlines be everywhere continuous, the perturba- 
tion stream function 7 must be continuous across the leading 
Mach wave (and hence must v ani sh as this Mach wave is approached 
from downstream) . This fact is obvious when equations (5) and (6) 
are considered to be integrated for UT, the path of integration 
crossing the Mach cone. 


1 ? 

2. On the body the boundary condition is that = — u^ R , 

where E is the radial distance of the leading edge from the axis 
of symmetry. Then, 


¥ 




( 8 ) 


where r-g(x) defines the body. 


These boundary conditions, together with equation (7), pro- 
vide a boundary -value problem of the first kind. 

Stream function for supersonic source . - The function 

satisfies differential equation (7). This function 
is the stream function for the well-known supersonic source. A 
distribution of such sources along the x-axis yields 
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¥ 


f(0) 


pz-pr 


'N 




u 


o \Jx-i) z -£' d r 2 
px-pr 

f*(|) V( z -0 2 -P 2 r 2 cL| 




U' 


J 


0 ) 


Substitution of equation (9) into equations (5) and (6) pro- 
vides expressions for u and v identical to those provided by • 
the distribution of source potentials (reference 2) . Substitution 
of equation (9) into the boundary condition given by equation (8) 
provides an integral equation for f(x). If a slender pointed 
body is under consideration, this integral equation (in •which B 
is set equal to zero) may be differentiated with respect to x 
and then, provided rg(x) is vanishingly small, may be solved to 
yield the relation 


f(x) 


i ^d_ r 

“ 2 U ® dx l_ 3 


r B ( x ) 


(x >0) 


= 0 j (x< 0) 


( 10 ) 


where the origin of coordinates is placed at the nose of the body. 
This relation was first used by von Karman (reference 2). 


Example; right circular cone, at zero yaw, semi vertex 
angle tan"^ a . - By use of the method of sources, equations (9) 
and (10) may be combined to yield the stream function for this flow. 
Alternatively, the solution may be obtained by direct integration 
of the differential equation, as follows: The solution may be 

written in the form ¥ = r^h(^), where the origin is taken at the 
nose and pr/x. This form for ¥ is chosen in order to insure 

that the flow will be "conical." Substitution in equation (7) 
yields an ordinary differential equation for h(r)). This equation 
may be solved to yield 
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vlier© flnri Cg are constants of integration to be determined 

from the boundary conditions. Inasmuch as the perturbation stream 
function must vanish on the bov wave (tj = l), C 2 must be zero. 
The bo und ary condition given by equation (8) provides that 



An identical result can' be obtained by use of the source - 
distribution technique. 


( 11 ) 


An interesting feature of the solution is noted: 

The exact formulation of boundary conditions yields the value 
of C;l given by equation (ll). 3he axial perturbation velocity u 
may be obtained from, the stream function as follows; 



p 1 

» - u CL C cosh x — 

00 

The corresponding potential problem may be solved in a similar 
manner, in which case the use of exact boundary • conditions pro- 
vides that 
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Although the foregoing linearized results are identical, the exact 
solutions differ, thus emphasizing the incorrectness, in general, 
of using exact boundary conditions in the small -perturbation theory. 

Use of stream function . - The use of the stream function has 
the primary advantage of providing a ready means of computing 
stream l ines away from the body; for example, in the design of a 
tunnel-wall insert to eliminate wall 'interference. Equations (2) 
and (4), and the assumption that the stream function is continu- 
ous across the leading Mach cone, may be used to yield the follow- 
ing equation of a streamline: 

¥ t “ | Ve 2 


or 


1 

2 


u r 

oo 


+ ¥ = - u r_, 
2 co E 


2 


where the streamline is identified by the radial coordinate r-g 
with which it enters the Mach cone. Eor example, the following 
equation is obtained for a streamline of the cone flow treated in 
the preceding section: 



If, in the course of calculations, it is necessary to find the 
family of streamlines associated with flow over a slender axially 
symmetric body at zero yaw, it would be simpler to formulate and 
to solve the problem in terms of the. stream function. 

A method of linearized characteristics based on equation (7) 
should yield the streamline family very quickly. Velocities could 
then be computed by forward differences . 


Treatment of Open-Uosed Bodies of Revolution 

An alternative numerical technique for open-nosed bodies, based 
on the use of the stream function produced by a source distribution, 
is developed in the following analysis: 
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The family of source distributions that gives polyn omi a l 
streamline slopes near the lip of the body is sought. An illus- 
trative open-nosed body is shovn in figure 1. The definition 

r B (x) = B + 8(x) 


where 8 << B, 

is introduced into equation (8), yielding the approximate relation 


¥(x,r B (x)) s - uES(x) 

Equating expressions (9) and (12) yields 

px-prg 


( 12 ) 


u^ESCx) 


U 


(x-£) f(£) [(x-£) 2 -p 2 r/] l/2 d£ (13) 


Equation (13) may be expressed in terms of the quantity p = x-prg, 
which is taken to be very small (thus restricting consideration to 
the vicinity of the lip) . The following approximate Equation 
results : 




M/^r * 


U 


o 


This is Abel's integral equation, the solution of which is, 
according to reference 6 (par. 11.8): 


nx 


f(x) = IkiM A. 


* y P dxj 


■ iiliS 

o yx- £ 


(14) 


whence, 

1. If 8(p)co p, then equation (14) provides that f(x) oo x 

2. If 8(p) CO p 2 , then f(x) CO x^/ 2 . 


1/2 


3. If 8(p)co p , then f(x) CO x' 


5/2 


and so forth. 
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A similar analysis will show that if f (x) co x (as in the 

# 3 /2 

method of reference 1), then 5(p) CD |i . Thus, as stated in the 
INTRODUCTION, the open -nosed tody associated with a linear source 
distribution has at the lip (jjl = 0) a vanishing slope and an 
infinite curvature. 


In view of the foregoing results, a sequence of solutions 
for ¥ are sought, each solution (1^) corresponding to a source 

distribution on the axis of symmetry, of the form x n “^/2 , -where 
n equals 1, 2, 3, ... 


These solutions can be combined to represent a given body, 
with the assurance that this representation of the body will be 
regular (that is, will have a Taylor expansion) near the lip. 


The following dimensiohless quantities are defined: 
ts ^ 


X 

_ r 
s= — 

7) = 

s 

Rp 

R 


t 


where 

g = 

¥ 

u B' 
00 







(15) 


and so forth 


By ^comparison of these definitions with equations (4), (5), and 


Sq> 



+ g 


(16) 


_u_ _ _1_ 1 5g 

^oo s 


at 2 ° s \ s 

Using the usual linearized form for the pressure coefficient yields 





2.1^ 
p2 s 5s 


(17) 


Near the lip of an open -nosed body, the two-dimensional char- 
acter of the flow insures the validity and the adequacy of equa- 
tion ( 17 ) . Far downstream of the lip, provided there are no 
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further streamline corners , the flow will approach that for the 
corresponding pointed -nose "body of revolution. Higher order terms 
appearing in Bernoulli's equation would then he of sufficient order 
of ma gni tude for inclusion in the formula for pressure coefficient, 
although a strict mathematical “basis for their inclusion would he 
lacking. Because this report places emphasis on the forward por- 
tion of the “body, additional discussion of this question would not 
he pertinent. 


In the notation cf equation (15), a sequence of stream, func- 
tions g n are to he obtained, each member of the sequence to cor- 
respond to a certain axial source-distribution function ^n* Thus, 
when definitions (15) are substituted into equation (9) and the 
variable of integration is suitably redefined, the following 
expression (aside from a constant) is obtained: 


rt-s 


Sn 




(18) 


be 


In accordance with the previous discussion. 


is chosen to 


f n (t) a lfct 


n-l/2 


(19) 


where n is an integer, and k n is an undetermined constant. 
Substitution of- equation (19) into equation (18) and integration 
according to the method set forth in appendix A provide the fol- 
lowing results for g n and, by superposition, a terminating series 
of N terms for g: 


n+l/2 

&a = V 5 BnW 


N 


,n + l/2. 


g . £ tt—'XW 


( 20 ) 


By use of equation (17) p 2 C may be expressed in terms of g^ 

P n 


*n 


( 21 ) 


J 
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where 


H^tj) S -2 \jl^ | [-T|K(o) + E(oj] 

S SS [ 2,)2 “ft* (K ‘ E) + Vi + "X-2 H n -3* ' 


T (ri) = 2 K(0) 

1 V 1 ^ 






-3- (4+5t|)K-5t]E 

[1+n 


n > 1 




T n-1 + ^ T n-2 + V3 + * ‘ • +T 1> + 
4(H n-2 +E n-3 + * * * +E l\ > n> 1 


J 


and 


K(O) H 


E(a) = 


-vr/2 

Jo 

prt/2 


dcp (l-o 2 sin 2 cp) 


-1/2 ^ 


dcp ( l-o 2 sin 2 cp) 


1/2 


-VS 






"V] 


( 22 ) 


(23) 


(24) 


(The symbols K and E are complete elliptic integrals of the 
first and second kinds, respectively . ) The H^q) anti T n (q) are 
tabulated in table I, for n = 1, 2, 3, 4, and 5. 


The equation for determining the constants k n is developed 
as follows : On the body, the boundary condition given by equa- 

tion (8) requires that 
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g 


2 - a B 


) 


Thus, equation (20) provides that 





vhere 


a B = 1 + c(t) 

1/2 1/2 

Because k-j_ corresponds to f (t) CD t and t is the 
only distribution function providing a nonvanishing streamline 
slope at the lip, k^ is determined solely by the initial body 
slope. Near q «= 1, let q = 1 - D, vhere V<<1, and thus 
tg s 1 + u. Equation (25) then becomes 

€ b -k-jH-L (1 - V ) (26) 

From equation (2 i). 


0«/2 



Thus, 

3 

-qK + E * q rtU 
or 

%(1- p) « - “= v 
2 V2 

By use of the fact that t « 1 + P on the body, equation (26) may 
be differentiated to provide 
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The remaining values af are computed as follows : 

1. Equation (25) is formulated at 17 -1 points on the "body 
(exclusive of the lip, because all the H n vanish for i) c 1). 

2. The resulting set of 17-1 linear algebraic equations is 
solved for k 2 , k 3 , . . . , kjj. 

2 

■When the constant has been determined, p C p is computed 

directly from equation (21) . A more detailed outline of numerical 
procedure will be subsequently presented. 

Linearized treatment of attached bow shock . - The first 
approximation to the strength of the bow shock for an open-nosed 
body may easily be obtained from g]_ because, in the linear the- 
ory, character is tics of the same family never intersect, thus pro- 
viding that the shock strength is a function only of the initial 

O 

body slope. Defining the shock strength as A(p C ) across the 
shock yields 

Q -l/2 

A(p 2 G ) = fcjt lim T x (ti) 

* 11— >1 


where q approaches 1 from below. 

From equations (23) and (24), 

pit/2 

dcp = JL 

V2 


lim T x (q) 
q -^1 


Vi” 


and applying equation (27), 
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or 


A 

9 - 1/2 - 1/2 
A(p^C p ) = 2 6«(l)t E = 2e* (l)s E 


4 ( 0 p ) 


- S' (PH) 
P 





(28) 


Equation (28) gives tie linearized rule for tie attenuation of tie 
slock strength, with increasing distance from tie axis of symmetry. 
(Note tlat when rg/R = 1, equation (28) gives tie two-dimensional 
result . ) 


In view of tie foregoing result, it slould te possille to pro- 
pose a formula for tie slape of tie low wave. ,Sud a formula 
slould, of course, le decked against a solution ly tie metlod of 
clarac ter is t ics in order to estatlisl its value. At first it would 
seem tlat a first-order determination of tie slape of tie low slock 
wave would le in conflict witl one of tie lasic results of tie 
first-order (linearized) tleory, namely, that all slock and expan- 
sion waves appear to le located at Mad surfaces . 


Tie exact ollique -shock-wave relations may, lowever, le used 
to slow tlat a first-order flow deflection tlrough a slock wave is 
associated witl a first-order variation of slock position from tlat 
of tie Mad surface. Hence, if formulas are available- for tie flow 
deflection and the pressure rise across tie slock, loti presumed 
accurate to first order, and if tley are consistent (to tie '"first 
order) witl tie requirements of tie ollique -slock relations, these 
pressure and deflection functions may le used to determine, to the 
first order, tie only remaining quantity involved in the slock 
relations, namely, the slock inclination. 

1. In order to slow tlat tie pressure -rise and flow -deflection 
formulas stand in tie proper relation, equation (28) is first 
modified to 

A “ &'(0B)^2Lj ^ (29) 

Equation (16) may le written as follows when equations (20) and 
(27) and a definition of sj as tie radial coordinate of tie 
streamline entering tie Mad surface are used: 
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1 2 2VT - , 

— s + — 1 — e* 

2 jr 


d)t 3/2 Hl ( n ) 


(Note that inasmuch as this is the equation of a streamline, s is 
to "be considered a function of t.) 


The derivative of the preceding equation may he evaluated at 
Tj = 1 to provide the flow deflection through the shock. The fol- 
lowing result is obtained: 


or 



e*(l)s 


- 1/2 




E 






(30) 


This result was previously obtained in reference 3, wherein an 
entirely different approach was used. Adding equations (29) and 
(30) yields 



(31) 


Equation (31) must now be shown to be consistent with the 
oblique -Bhock relations. The angle between the shock wave the 
axis of symmetry, in the original (x,r) coordinates, is defined 
to be 0. This definition may be introduced into equation 4.22 
of reference 7, which then may be written as 


tan 0 = - ^ = - 


Au 


(vsu) (tL 


or. Inasmuch as tan 0 « 1/3, 


A #- 

1 u 

is -®- 


21 


q=l 
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which, clearly reduces to the first-order relation given in equa- 
tion (31). 

2. The shock-wave inclination may now he determined.' Equa- 
tion 4.29 of reference 7 may he written as 

\ 

1 + I M 2 A(Oj = — M 2 sin 2 6 - £li 
2 CO ' P' 00 -V+l 


1 + ^ M„ 2 4(C p ) 
tan 9 «= —g — 

P A(C p ) 


From equation (28), 


A(O p ) - | S' (SR)U 


and thus, because A(C^) is a sma ll quantity. 


. tan 0 s i + Zli 5J2L A(Cj 


P 8 p 3 


, q 1 7+1 . ^oo 

9 = + Z__ 

P 


6’ (PH) 


m 

\E ) 


If 0jj is defined at the shock inclination (fig. l) in the coordi- 
nates introduced in equations (15), equation (32) becomes 

4 

7+1 M oo . . -l/2 

tan 8^ s 1 + 1 ^ €* (l)sjj (33) 

Equation (33) is exact' at an infinite distance from the body. 
The error at the lip of the body should indicate the magnitude of 
error at finite distances from the surface. For a free -stream Mach 
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number of 1.50 in air ( y = 1.4) and an initial "body slope of 
5 percent, tlie exact oblique -shock relations yield a value of 
tan 9 «= 1 at the lip; whereas equation (33) gives tan 9 = 0.988. 
These results are to he compared with a Mach angle equal to 
tan"-*- 0.894. At the lip, equation (33) is therefore in error by 
about 1 percent in predicting the shock angle and by about 12 per- 
cent in predicting the difference between shock and Mach angles. 


METHOD FOE ERESSDBE CALCULATION AHD EXAMPLE 

The functions T n and H n , appearing in equations (21) and 
(25), respectively, are tabulated in table I. The reciprocal 
sequence of tj provides approximately equal intervals of tg(Tj) 
for a slender body, thus making it unnecessary to perform inter- 
polations. The following numerical procedure is recommended: 

1. Select the number (H-l) and location (t^) of the points 
(exclusive of the lip) through which the streamlines will be 
required to pass. Define these points by values of tj for which 
the functions in table I are given. 

2. Compute tg(T^) from the equation of the body and tabulate 

Vn)' 172 , Mm) 1/2 , • • • . t B (ii) IM/2 - 

3. Compute kj_ from the initial slope of the body, using 
equation (27). 

4. Formulate equation (25) at the H-l "matching points," 
and solve the IT-1 simultaneous equations for k g , k 3 , . . . , k^. 

2 

5. Use the k n thus determined to compute ^ C p at any- 
desired location, using equation (21). 

This procedure has been carried out for a truncated cone of 
slope e'(l) = 0.0523, matching streamlines at tj = l/2, l/3, 
and l/4, and at tj = l/z and l/3. The results appear in figure 2. 
Comparison is made with the result obtained by the method of refer- 
ence 1 for a truncated cone of the same slope, using approximately 
eight integration points in the range of t shown (reference 8). 

An additional curve shown in figure 2 was obtained by applying the 
boundary condition on streamline slope at tj n l/2, l/3, and l/4 
on the body, and computing C-p, as was previously done. Appen- 
dix B contai n s the pertinent analysis and procedure, which gives 
rise to the function ) plotted in table I. 
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In order to determine the accuracy attained hy using a certain 
set of matching points, equation (25) may easily he inverted and 
solved for € at intermediate points. These values of c should, 
of course, fall on the surface of the given body, -with an error of 
the order of e 2 . This inversion may he carried out with suf- 
ficient accuracy hy assuming that tjj(q) at the intermediate 
points is given hy the equation of the body. 

As a guide to the selection of matching points, the following 
generalizations are offered: 

1. The farther downstream the computations are to he extended, 
the denser must he the matching -point distribution. (In the 
method of reference 1, the density of integration stations must he 
greatest near the lip.) This question of convergence. will he sub- 
sequently discussed. 

2. Very few points are needed near the lip. 

3. Pressure values cannot he computed with adequate accuracy 
at points downstream of the last matching point. 


DISCDSSIOn OF METHOD FOR PRESSURE CALCULATION 

The use of the half -power source distribution on the axis of 
symmetry provides a family of stream functions that are easily 
related to the forward portion of a given body. The component 
stream functions of equation (25), however, became infinite at 
y «= 0 and, as computations are extended downstream, satisfactory 
convergence is obtained with increasing difficulty; that is, the 
matching -point density must he greater. Hear the lip of an open- 
nosed body, the present method therefore converges rapidly; whereas 
the method of reference 1 converges slowly (the computations of 
reference 1 do not automatically provide the two-dimensional solu- 
tion at the lip, hut rather converge to that solution as more and 
more integration stations are used) . Toward the aft portion of 
the body, the reverse is true. 

Good agreement between the method proposed herein and the 
method of reference 1, on the forward portion of the body, is shown 
In figure 2. Farther downstream, the result begins to diverge from 
the proper one - for the reason previously mentioned. In order to 
extend the solution downstream, more matching points would have to 
be provided. 


u • 
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Of course, the stream function is not essential to the present 
technique, nor is it the cause of the downstream divergence of the 
method. For the case computed, the curve in figure 2 obtained by 
matching streamline slopes (as is customary in analysis based on 
the velocity potential) diverges from the proper result farther 
downstream than does the corresponding curve obtained by matching 
streamline locations. 


CONCLUSIONS 

Stokes' stream function for supersonic linearized flow may be 
used to advantage whenever it is desirable to deal directly with 
streamlines . 

Use of half -power source distributions on the axis of symmetry 
gives a family of stream functions having streamlines regular near 
the lip of an open-nosed body of revolution. These streamlines, 
being regular, are easily matched to the forward portion of a given 
body. 


First-order relations have been obtained for the attenuation 
of shock strength and flow deflection outward along the shock and 
for the corresponding shock-wave inclination. These results should 
be checked against experiment or the method of characteristics 
before being applied to practical problems. It is felt that these 
relations would probably conform rather well to physical fact, pro- 
vided the entrance slope Is small. These relations contain an 
error of the order of [&* (pE)J 2 at the lip of the body (where 5 
is a small quantity depending on the body shape, p is the cotan- 
gent of the Mach angle, and E Is the radial coordinate of the 
lip of the bojLy), and become exact as the radial coordinate r 
approaches infinity. It thus seems plausible that the error in 
the formulas for shock strength and inclination would be of second 
order over the entire shock front. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, December 5, 1949. 
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APEEKDU A 

METHOD OF EVALUATING g^ 

Equation (18) may "be written, after the substitution of - equa- 
tion (19), as 


nt-s 




°i)j 0 


[(t-T) 


2 211/2 n-3/2 

-S T 


4T (Al) 


The following change of variable of integration is introduced: 

-(i^f ' 

This substitution yields elliptic integrals in Legendre * b canonical 
forms. From reference 6 (par. 22.72), 


r\i 


UO 


ay jii-y 2 )(l-ay 2 )J ^ = K(a) 


(A2) 


ni 


1/2 


y 2 dy Jj(l-y 2 ) (l-o 2 y^)j . = i [K( a )-E(a)] (A3) 


UO L -* <r 

differentiating and Integrating the expression 


2n-l T/i 2w n 2 27] 
y ]_(i-y )(i-o y )| 


1/2 


yields the following recursion equation: 
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ni 


r\i 


(2n-l) 


y 2n “ 2 dy j( 1-y 2 ) ( l-^T 2 ! y 2 )| ^ -2n( l+c ) ^ 

[. . *r 


U< 


2n, 
y dy 


E-3 


- 1/2 


(2n-l)o 


2n+2 
7 dy 


(M) 


Equations (A2), (A3), and' (A4) permit equation (Al) to "be expressed 
in terms of K(a) and E(cr). 


In computing velocities, use is found for the following equa- 
tions (reference 6, par. 22.736): 


0 


dE(c) 



[jS(cr)-K(a)] 


dg(q) 

do 2 


— r ^TE(a)-K(a) 

2 o 2 _l-a 2 
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APPENDIX B 


BOUNDARY CONDITION ON STREAMLINE SLOPE 

In the numerical method developed for the determination of 
pressure on an- open-nosed body, a Boundary condition (equation (8)) 
on the stream function has Been applied. Alternatively, a Boundary 
condition on streamline slope may Be used as follows: The Boundary 

condition is 



6'(x) 


vhich Becomes, after definitions (5) and. (6) are introduced. 


. _ ZL^V 

s B \dt t drj 


€'(t) 

0 


(Bl) 


Equations (15), (17), and (21) enaBle to Be expressed in terms 

dij 

of T n (i}), and equation (20) provides a relation Between 
and H n (Ti). Thus, equation (Bl) Becomes 


dt 


and 


‘■(tsW)- i 1* t B (T)) M/2 Q tt (t|) 


(B2) 


where 


Values of the function (^(tj) may Be found in taBle I. 

The method for pressure calculation using the Boundary con- 
dition on streamline slope is precisely the same as that previously 
described, with the exception that equation (B2) has the role here- 
tofore played By equation (25). 
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TABLE I - FUNCTIONS USED FOR COMPUTING STREAMLINES AND PRESSURES 


1/1 

-Hi 

-h 2 

-h 3 

-h 4 

-h 6 

1 

0 

0 

0 

0 

-0 

1.25 

.20766 

.03054 

.00475 

.00032 

-.00080 

1.5 

.32871 

.07922 

.02042 

.00392 

-.00194 

2 

.46001 

.16453 

.06697 

. 02903 

-.01318 

3 

.56426 

.26112 

.13912 

.07914 

-.04681 

4 

.60436 

.30799 

.18195 

.11505 

-.07560 

5 

.62480 

.33451 

.20873 

.13987 

-.09763 

6 

.63624 

.35063 

.22603 

.15679 

-.11344 

7 

.64346 

.36132 

.23797 

.16893 

-.12519 

8 

.64849 

.36888 

.24667 

.17802 

-.13426 

9 

.65204 

.37437 

.25311 

.18487 

-.14121 

10 

.65454 

.37845 

.25801 

.19019 

-.14669 

12 

.65764 

.38390 

.26475 

. 19770 

-.15460 

14 

.65986 

.38758 

.26935 

.20286 

-.16011 

16 

.66114 

.39001 

.27246 

.20642 

-.16400 

i/u 

Tl 

t 2 

t 3 

T 4 i 

?5 


1 

1.25 

1.5 

2 

5 

4 

5 

6 

7 

8 
9 

10 

12 

14 

16 


2.22171 


1.25 

2.41109 

. 73505 

.18679 

.04786 

.01895 

1.5 

2.57167 

1.32770 

.56925 

.24245 

. 12785 

2 

2.83099 

2.22958 

1.41591 

. 83058 

.46526 

3 

3.21120 

3.48473 

, 2.98658 

2.35334 

1.77371 

4 

3.48788 

4.37032 

4.26209 

3.81094 

3.26231 

5 

3.70498 

5.04664 

5.28699 

5.06011 

4.62728 

6 

3.88312 

5.59980 

6.15194 

6.15842 

5.88378 

7 

4.03528 

6.06782 

6.89588 

7. 12501 

7.01890 

. 8 

4.16621 

6.46680 

7.53546 

7.96738 

8.02372 

9 

4.28338 

6.82305 

8.11157 

8.73563 

8.95364 

10 

4.38734 

7.13981 

8.62701 

9.42871 

9.80062 

12 

4.56878 

7.69240 

9.53250 

10.65814 

11.32064 

14 

4.72204 

8.15524 

10.29267 

11.69685 

12.61471 

16 

4.85539 

8.55928 

10.95993 

12.61379 

13.76563 

l/l 

<*L 

^2 

Q3 

Q4 

$5 


1.11086 

1.35380 

1.59683 

2.08778 

3.07436 

4.06215 

5.05650 

6.04978 

7.04457 

8.04227 

9.04050 

10.03747 

12.02788 

14.02571 

16.01909 


•0 

.38946 

.73965 

1.38005 

2.53920 

3.62619 

4.68604 

5.72611 

6.75652 

7.78178 

8.80238 

9.81824 

11.83750 

13.85656 

15.86788 


0 

.09550 

.29696 

.82277 

1.95852 

3.08006 

4.18148 

5.25930 

6.32284 

7.37773 

8.42361 

9.46170 

11.51667 

13.56575 

15.60026 


0 

.02094 

.10728 

.46892 

1.46061 

2.54727 

3.65309 

4.74654 

5.83024 

6.90668 

7.97255 

9.02999 

11.11987 

13.19793 

15.25642 


0 

. 00208 
.02661 
-.02866 
-.47675 
-1.25541 
- 2.22210 
-3.25320 
-4.31846 
-5.40596 
-6.49248 
-7.57792 
-9.73193 
-11.87794 
-14.00182 



Axis of symmetry 



Figure 1. - Illustrative open-nosed body. 










